home *** CD-ROM | disk | FTP | other *** search
/ Aminet 2 / Aminet AMIGA CDROM (1994)(Walnut Creek)[Feb 1994][W.O. 44790-1].iso / Aminet / util / misc / Fudgit233.lha / Source / examples / bessel.ft < prev    next >
Encoding:
Text File  |  1993-12-14  |  626 b   |  33 lines
  1. # This file plots first bessel functions using GNUPLOT
  2. pmode
  3.     set data style line
  4.     set nokey
  5. fmode
  6. set data 500
  7. cmode
  8.     x=0
  9.     X=x++
  10.     X*=(20/data)
  11.     J1 = besjn(1, X)
  12.     J2 = besjn(2, X)
  13.     J3 = besjn(3, X)
  14.     J4 = besjn(4, X)
  15.     J5 = besjn(5, X)
  16.     X2=X+3
  17.     Y1 = besyn(1, X2)
  18.     Y2 = besyn(2, X2)
  19.     Y3 = besyn(3, X2)
  20.     Y4 = besyn(4, X2)
  21.     Y5 = besyn(5, X2)
  22. fmode
  23. save vec X J1 J2 J3 J4 J5 $Tmp.1
  24. save vec X2 Y1 Y2 Y3 Y4 Y5 $Tmp.2
  25. pmode
  26.     plot '$Tmp.1' us 1:2, '$Tmp.1'  us 1:3, '$Tmp.1' us 1:4,\
  27.     '$Tmp.1' us 1:5, '$Tmp.1' us 1:6
  28.     pause 3
  29.     plot '$Tmp.2' us 1:2, '$Tmp.2'  us 1:3, '$Tmp.2' us 1:4,\
  30.     '$Tmp.2' us 1:5, '$Tmp.2' us 1:6
  31.     pause 3
  32. fmode
  33.